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Abstract 



Using the quark mass density- and temperature dependent model, we have 
studied the thermodynamical properties and the stability of strangelet at 
finite temperature. The temperature, charge and strangeness dependences on 
the stability of strangelet are investigated. We find that the stable strangelets 
are only occured in the high strangeness and high negative charge region. 

PACS number: 12.39.Ki, 21.65.+f, ll.10.Wx, 25.75.-q 
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I. INTRODUCTION 



The study of small lumps of strange quark matter, called strangelets, plays an important 
role for research the quark gluon plasma (QGP) in recent relativistic heavy ion collision 
(RHIC) experiments. The reason is that although many signatures of QGP such as J/^J 
suppression, strangeness enhancement, thermal dilepton electro magnetic radiation, etc., 
have been found [|I],[|, but it is still ambiguous because these signatures can also be ex- 
plained by hadron gas ||. To search an unambiguous signature of QGP is the key for RHIC 
experiments. The strangelet, as was argued by Greiner et. al f|, is a good candidate which 
could serve as an unambiguous signature for the QGP. 

The essential problem for detectability of strangelet in RHIC experiments is to study its 
stability during the formation of QGP. Employing MIT bag model, many authors discussed 
this problem They came to the same conclusion that the electric charge and the 

strangeness fraction are of vital important to the experimental searches of strangelets because 
these two elements affect a change of stability remarkably. But the effect of charge on 
stability is different for different authors. Jaffe and his co-workers [|5|-[7| argued that the 
strangelet is slightly positive charged. They considered strong and weak decay by nucleon 
and hyperon emission together and concluded that the stable strangelets will have a low but 
positive change to mass ratio. Contrary to Jaffe et. al, after considering the initial condition 
of possible strangelet production in RHIC carefully, Greiner and his co-workers || argued 
that strangelets are most likely highly negative charged. Here we hope to emphasize that 
above discussions are limited in the framework of MIT bag model and at zero temperature. 

Since the quark deconfmement phase transition can occur at high temperature or/and 
high density only, it is of interest to extend the investigation for the stability of strangelets 
to finite temperature. This is the objective of this paper. 

However, MIT bag model is a permanent quark confinement model because the confined 
boundary condition does not change with temperature. In principle, one can not use this 
model to study the phase transition of QCD directly. The best that we can do is to employ 
a model which can almost reproduce the properties of strange quark matter obtained by 
MIT bag model, but in which the quark confinement is not permanent. The quark mass 
density- and temperature-dependent model (QMDTD) [§,[K| suggested by us is one of such 
candidates. 

Many years ago, a quark mass density- dependent model (QMDD) was suggested by 
Fowler, Raha and Weiner |TT|] and then it was employed by many authors to discuss the 
properties of strange quark matter fl2|-|T5||. According to the QMDD model, the masses of 
u, d quarks and strange quarks (and the corresponding anti-quarks) are given by 

m q = - — , (q = u,d,u } d), (1) 
3n B 

m StS = m s0 + - — , (2) 
3n B 

where % is the baryon number density, m s0 is the current mass of the strange quark and B 
is the vacuum energy density. As was proved by ref. [13], the properties of strange matter in 
the QMDD model are nearly the same as those obtained in the MIT bag model. In fact, it 
is not surprised if one notices that the confinement mechanism of MIT bag model is almost 
the same as that of QMDD model ||. 
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But when we employ the QMDD model to discuss the properties of strangelet, a lot 
of difficulties emerge At first, the radius of strangelet decreases as the temperature 

increases; Secondly, it can not mimic the correct phase diagram of QCD because the tem- 
perature T tends to infinite when ub — > 0. To overcome these difficulties, we suggest a 
QMDTD model 0,0. Instead of a constant B in QMDD model, we argue that B is a 
function of temperature and introduced an ansatz JTO] 



B(T) = B 



,0<T<T c , (3) 



B(T) = 0,T>T c , (4) 

where B is the vacuum energy density inside the bag (bag constant) at zero temperature, 
T c = 170MeV is the critical temperature of quark deconfinement phase transition, and a, b 
are two adjust parameters. Since B is zero when T = T c , a condition 

l-a + b = (5) 

is imposed and only one parameter a can be adjusted. As pointed out by ref. ||10|| , in order 
to satisfy two physical conditions of strangelet at finite temperature, namely, (1), the radius 
of the strangelet must increases when the temperature rises, (2), the energy of the strangelet 
must increases when the temperature rises, the parameter a is restricted in a small range: 

0.65 < a < 0.8. (6) 

In this paper, we fix the value of two parameters a, b in this suitable range 

a = 0.65,6= -0.35. (7) 

With these parameters set, we use the QMDTD model to discuss the thermodynamical 
properties of strangelet, especially, to investigate the stability of strangelet via strong hadron 
emission and weak hadronic decay. We will study the effects of temperature, charge and 
strangeness fraction on the stability and hope that our study can have impact to the de- 
tect ability of the strangelet. 

The organization of this paper is as follows. In the following section, we give the formu- 
lae of thermodynamical calculations for QMDTD model. The results of thermodynamical 
properties of strangelet are presented in section 3. In section 4, we will discuss the stability 
of strangelet via possible strong and weak decays. The last section is a summary. 

II. THERMODYNAMICAL FORMULAE 

To calculate the dynamical and thermodynamical quantities of the strangelet, we must 
look for the density of states first. The density of states of a spherical cavity in which the 
free particles be contained can be expressed as 

m = ( 8) 
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where N(k) is the total number of particle states and can be written in terms of dimensionless 
variable kR as 



N(k) = A{kRf + B{kRf + C(kR), 



(9) 



where R is the radius of the bag. The three terms of the right hand side of Eq.@ refer 
to the contributions of the volume, surface and curvature, respectively. The coefficients A, 
B and C are expected to be very slow varying functions of kR, and their expressions are 
model dependent. For MIT bag model, these coefficients have been obtained by numerical 
calculations in our previous paper [|16j. The volume term A is a constant, and it has the 
value of 



A 



2l 

9tt ! 



(10) 



where g is the total degeneracies. For example, it is the total number of spin and color 
degrees of freedom for a quark with flavor treated separately. The surface term B is 



B\-\=tL 



k 



2tt 



1 -'t)' 



tan"' - - - -1 



m 



n 



k) 



(11) 



where m is the mass of quark. Equations ( |T0"D and (|TT|) are in good agreement with the ones 
given by multi- reflection theory |T^JT3] • The curvature term C can not be evaluated by this 



theory except for the two limiting case: m — > and m — > oo. Madsen proposed that |]T9 

k 



~k' 



9_ 

2tt 
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rn 



ilk 
2m 



(12) 



But as was pointed out by ref. [IB], the beat fit of numerical data for curvature term is 

1.45 



J' 
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3.42 ( 



m 



6.5 +100 



(13) 



Now we are in the position to calculate the thermodynamical quantities of strangelets 
for QMDTD model. The thermodynamical potential Q is 



n = E"* = -E 



(2tt) 3 



dk 



dk 



In 1 



-P(ei(k)—m) 



)■ 



dNi 
dk 



(14) 



is the 



where i stands for u, d, s (or u,d,s ) quarks, = 6 for quarks and antiquarks. 

density of states for various flavor quarks, it is given by eqs.(p[)-(P^). jii is the corresponding 
chemical potential (for antiparticle m = — /ij). 



Si(k) = \J mf + k 2 



(15) 



is the single particle energy and mj is mass for quarks and antiquarks. 
According to the QMDTD model [pyTOf, the masses of quarks are 
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B, 



rn 



o 



u'u.d.d 



3n B 



T 



a — )+b — 



T 



T, 



, < T < T r 



m u ,u,d,d = 0, T > T, 



CI 

3n B 



T I 



, < T < T r 



0, T > T r 



(16) 



(17) 



where m sQ is the current mass of the strange quark matter, % is the baryon number density 

n B = A/V, (18) 
4 

and A is the baryon number of the strangelet, V = g 71 "^ 3 is the volume of the strangelet. 
Using the standard statistical treatment, and noticing that Q is not only a function of 
temperature, volume and chemical potential, but also of density, it can be proved that the 
total pressure p and the total energy density e are given by HHiHi 



V = - 



1 d(Q/n B ) 



V d(l/n B ) 



-T7 + 



V V dn B 



(19) 



Q ^ T dQ 

£ = h) UiUi — — 

v i V dT 

The number density of each particle can be obtained by means of 

1 dtt 



(20) 



Hi = - 



V d/j>i 



(21) 



T,n B 



At finite temperature, we must include the contributions of the anti-particles, therefore, the 
baryon number for i quark is given by 



AN t = (m -n^xV = [°°dk^ ( -— - -— - 

V l! (2tt) 3 7o dk \exp\J3(ei-fH)] + l expf/% + ^)] + 1, 



The strangeness number S of the strangelet reads 

S = AN S , 

the baryon number A of the strangelet satisfies 



(22) 



(23) 



A=-(AN u + AN d + AN s ). 



The electric charge Z of the strangelet is 



Z = ^AN U - ^AN d - ^AN S . 



(24) 



(25) 



At finite temperature, the stability condition of strangelets for the radius reads 

where the free energy F of strangelet is, 

F = E-TS, (27) 

E — eV is the total energy, and 



(28) 



is the entropy 

Given strangeness number, baryon number and electric charge, for any strangelet, 
we can calculate chemical potentials for u, d, s quarks self-consistently by the equations 
(p3|).(p^),(p5D and (|26|). Then the thermodynamic potential, free energy and the stable 
radius of the strangelet at finite temperature can be obtained self-consistently 

III. THERMODYNAMIC AL PROPERTIES OF STRANGELET 

The numerical calculations have been done with the parameters set: 

B = i70MeVfnT 3 , m s0 = 150MeV, T c = 170MeV, (29) 

and the possible area for strangelet in our calculations is chosen as 

S > 0, (30) 
Z > -A, (31) 
S + Z<2A. (32) 

We calculate the free energy of the strangelet first. Figure 1 shows the free energy per 
baryon F/A in a S — Z plane for all possible strangelets with A = 5 and T = 20MeV. This 
is a set of downward protruding curves. The corresponding curves but for T = 50MeV are 
shown in figure 2. Comparing these two figures, we find that the positions of dots appeared 
in figure 2 are lower than the corresponding positions (with the same strangeness number S 
and charge Z) in figure 1. For example, for S = 3 and Z = —2, the free energy per baryon of 
the dot is F/A = 1064.9MeV in figure 1, but reads F/A = 898.9MeV in figure 2. For fixed S 
and Z, the free energy of the strangelet decreases when temperature increases. This result 
is reasonable because the entropy increases with the temperature, the term TS in eq. (p7|) 
increases considerably and F decreases. We will see in the next section this result affects on 
the decay of the strangelet remarkably. 

To illustrate our result transparently, we draw the F/A vs S curves for fixed A = 5, Z = 
—5 but different temperatures T = 20MeV and 50MeV in figure 3, respectively. We find 
that the minimum of these two curves are located at the same strangeness S — 4 but with 
different free energies.. The whole curve for T = 20MeV is located on the upper position 
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of the curve for T = 50MeV. We also draw the F/A vs S curve for A = 10, Z = —10 and 
T = 50MeV in figure 4. Compared with figure 3, the shape of the curve in figure 4 almost 
does not changes except the position of the minimum is changed to S = 8. These results 



are similar to that one given by |20 



Now we turn to study the charge and strangeness dependences of the radius of strangelet. 
We draw F/A vs R curves with fixed A — 5 ,Z — 1 and T = 50MeV but different S = 2, 6 
and 9 in figure 5, respectively. The stable radiuses given by Eq.(|2TS|) for different strangeness 
are different. Figure 5 shows that the stable radius changes from 1.64fm to 1.615fm when 
strangeness number changes from 2 to 9. The same curves for A = 5 ,Z = —4 and T = 
50MeV but different S — 2, 14 are shown in figure 6. We see that the stable radius become 
1.67fm for S = 2 and 1.62fm for S = 14. Therefore, we come to a conclusion that the stable 
radius of strangelet decreases when S increases. 

To illustrate the charge dependence of the stable radius, we draw F/A vs R curves with 
fixed A = 5 ,S = 2 and T = 50MeV but different electric charge Z — 1,5 and 8 in figure 
7, respectively. The stable radius increases from 1.64fm to 1.69fm when the electric charge 
rises from 1 to 8. The same curves for A = 5 ,S — 9 and T = 50MeV but different charge 
Z — 1, —4 are shown in figure 8. We see that the stable radius changes from 1.615fm 
to 1.605fm when Z decreases from 1 to —4. We find that the stable radius of strangelet 
increases with electric charge . 

IV. STABILITY OF THE STRANGELET 

In this section, we follow the line of ref. M to investigate the stability of the strangelet 
and extend their study to finite temperature by using the QMDTD model. 

A. Strong decay and unstable strangelet 



As was pointed out by refs. [pTp2| , small clusters of strange matter are most favoured 



for detection. As two examples, hereafter we study two cases with A = 5 and A = 10, 
respectively. 

Instead of the binding energy at zero temperature, we calculate the free energy per baryon 
of the possible strangelet at finite temperature first. Figure 9 and 10 show the free energy 
per baryon of all possible strangelets as a function of the strangeness S at the temperature 
T = 50MeV but for A = 5 and A = 10, respectively. The solid lines in these two figures 
connect the masses of the nucleon, A, 5 and Q. As a first cut for potential candidates of 
stable strangelets, these ones lying above this line can (or probably will) completely decay 
to the pure hadron state via strong processes, and only these ones beneath the line will be 
possible for metastable or stable strangelets ||. 

We consider the influence of temperature now. As shown in last section, the free energy 
per baryon increases when temperature decreases. The positions of the dots will raise and 
many dots will across the line and become unstable when temperature decreases. Figure 11 
shows this result clearly. The figure 11 is the same as figure 9 except for temperature T = 
20MeV. Comparing these two figures, we find many dots across the line when temperature 
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decreases from 50MeV to 20MeV. At low temperature, more strangelets can decay to the 
pure hadronic state by strong process. 

Here we must emphasize that although any strangelet initially formed under this line 
can not decay to a pure hadron state, it is still possibly unstable because it can decay to 
a hadron and another strangelet changing baryon number, strangeness number and charge 
||. We will look for possible strong decays, i.e. single baryon {n, p, A, S + , E~, 5°, H~ and 
fl} emission and mesonic decays at the finite size configuration at finite temperature. 

The baryon number A, strangeness number S and electric charge Z are conserved in the 
strong process. A general expression of a strong baryon decay for a strangelet Q(A, S, Z) 
can be written as 

Q(A, S, Z) -»• Q(A-1,3- S x , Z-Z x )+ x(l, S x , Z x ). (33) 

And this process is allowed if the free energy balance of the corresponding reaction is 

F(A, S, Z) > F(A — 1,S — S x , Z - Z x ) + m Xj (34) 

where F stands for the total free energy of the strangelet and x stands for a baryon with 
strangeness number S x and electric charge Z x . Our results of numerical calculation are 
shown in figures 12, 13 and 14. For A = 5 and T = 50MeV, all strangelets including the 
ones lying above the lines of figure 9 and the ones which satisfy the inequality Q54]) are 
drawn by the circles in figure 12. These strangelets can strong decay and are unstable. The 
unstable strangelets for A = 10 and T = 50MeV are drawn by circle in figure 13. Figure 
12 and figure 13 show that the strangelets with small strangeness number S and situated in 
the left side of the figure are unstable. 

To study the temperature effect, the same figure for A = 5 and T = 20MeV are shown 
in figure 14. Comparing figure 12 and figure 14, we find that many strangelets become 
unstable when temperature decreases. Lower temperature is in favour of the strong decay 
of the strangelet. 

B. Weak decay and metastable strangelet 

According to the definition of Greiner et. al 0, "a strangelet is called metastable in the 
following if its energy lies under the correspond (free) hadronic matter of the same baryon 
number, charge and strangeness, and if it can not emit a single hadron or multiple hadrons by 
strong processes." At finite temperature, instead of energy, we use free energy. A metastable 
strangelet can then only decay via weak process like the nonleptonic (hadronic) decays. In 
this sub-section, we study all possible weak hadronic decay for metastable strangelets. 

In the weak process, the baryon number A and the electric charge Z are conserved, but 
the strangeness number S is not conserved. For the weak decay, AS = ±1. Therefore, a 
general expression of a weak baryon decay for a strangelet Q(A, S, Z) can be written as 

Q(A, S, Z) -> Q(A - 1, S - S x - 1, Z - Z x ) + x(l, S x , Z x ). (35) 

This process is allowed if the free energy balance of the corresponding reaction is 

F(A, S, Z) > F(A — 1,S— S x -l,Z-Z x ) + m, (36) 



S 



The metastable strangelets which satisfy the inequality fl36|) are shown in figure 12, 13 and 
14 with filled circle. We find that they are almost situated in the area with negative charge 
(or slightly positive charge) and high strangeness number. In particular, high strangeness 
is in favour of the stability of strangelet. As shown in figures 12, 13 and 14, a minimal 
strangeness S c above which (S > S c ) strangelets are metastable or stable exists. The values 
of the minimal strangeness S c are 5 in figure 12 (A = 5,T = 50MeV) and 10 in figure 13 
(A = 10, T = 50MeV). Defining strangeness fraction f s as 

/. = j, (37) 

Greiner and his co-workers had predicted that there exists a critical f sc at zero tempera- 
ture and pointed that the metastable or stable strangelets could only be found above this 
value f s > fsc. Our results strongly support their prediction and extend it to the finite 
temperature. We obtained f sc — 1 at T — 50MeV. This value does not depend on the 
baryon numbers. But as shown in figure 14, f sc depends on the temperature. It becomes 

f sc = — = 2 at T = 20MeV. 
5 

Now we turn to discuss the charge dependence of the stability of strangelet. As shown 
in figures 12, 13 and 14, higher negative charge is in favour of the metastable and/or stable 
strangelets. A maximal charge, which is Z m = 1 in figure 12 (A = 5,T = 50MeV), Z m = 6 
in figure 13 (A = 10, T = 50MeV) and Z m = -3 in figure 13 (A = 5,T = 20MeV), 
respectively, exists. The strangelets would be stable or metastable when Z < Z m . Defining 
charge fraction f z as 

/, = f, (38) 

We find that the maximal charge fraction depends on not only the temperature but also 

1 6 

the baryon number A. For example, f zm equals to - for A = 5 and — for A = 10 when 

5 10 

T = 50MeV. 



C. Stable strangelet 

The strangelet which is stable against both strong and weak decay is called stable 
strangelet. The stable strangelets are shown in figure 12, 13 and 14 by filled squares. 
We find that only a few strangelets can not decay via strong and weak reactions and be 
completely stable. For example, as shown in figure 12, for A = 5 and T = 50MeV, in 
total 128 strangelets, only 4 strangelets are stable [Z = —3,5 = 12), {Z = —4,5 = 13), 
(Z = -4, S = 14) and (Z = -5,S= 15). 

Finally, we hope to emphasize that the high negative charge and the high strangeness 
are in favour of the stable strangelets obviously at finite temperature. The stable strangelets 
in figures 12, 13 and 14 are all highly negative charged. The conclusion given by Greiner et. 
al at zero temperature is still correct at finite temperature. 
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V. SUMMARY 



In summary, by using the QMDTD model, we have studied the thermodynamical prop- 
erties and the stability of strangelets at finite temperature. We obtain: 

(1) , For fixed strangeness and charge, free energy per baryon decreases as the temperature 
increases. The stable radius of the strangelet decreases when the strangeness increases or 
the charge decreases. 

(2) , The higher temperature is in favour of the stability of the strangelets. Comparing 
figure 12 with figure 14, we see clearly that the total area of metastable and stable strangelets 
expands when the temperature increases. 

(3) , The higher negative charge and higher strangeness number are in favour of the 
stability of the strangelets. The stable strangelets are highly negative charged and have 
high strangeness. 
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VII. FIGURE CAPTIONS 

Figure 1. The dots stand for the free energy per baryon F/A in a S — Z plane for all 
possible strangelets with A = 5 and T = 20MeV. 

Figure 2. The same as figure 1, but for the temperature T = 50MeV. 

Figure 3. The free energy per baryon F/A as functions of strangeness number S with 
A = 5 and Z = — 5, the two lines represent the temperature T = 20MeV and 50MeV, 
respectively. 

Figure 4. The free energy per baryon F/A as a function of strangeness number S with 
A = 10 , Z = -10 and T = 50MeV. 

Figure 5. The free energy per baryon F/A as functions of radius R with A = 5 , 
Z — 1 and T = 50MeV, the three lines represent the strangeness number ,5 = 2,6 and 9, 
respectively. 

Figure 6. The free energy per baryon F/A as functions of radius R with A = 5 , Z = —4 
and T = 50MeV, the two lines represent the strangeness number S — 2 and 14, respectively. 

Figure 7. The free energy per baryon F/A as functions of radius R with A = 5 , S = 2 
and T = 50MeV, the three lines represent the electric charge Z = 1, 5 and 8, respectively. 

Figure 8. The free energy per baryon F/A as functions of radius R with A = 5 , S = 9 
and T = 50MeV, the two lines represent the electric charge Z = —4 and 1, respectively. 

Figure 9. The dots stand for the free energy per baryon F/A with various strangeness 
numbers S for all possible strangelets with baryon number A = 5 at temperature T = 
50MeV. The masses of nucleon, A, S and fl are represented by the filled squares, respectively. 

Figure 10. The same as figure 9, but for baryon number A = 10. 

Figure 11. The same as figure 9, but for temperature T = 20MeV. 
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Figure 12. The electric charge Z as a function of the strangeness number S for unstable 
strangelets (open circles), metastable strangelets (filled circles) and stable strangelets (filled 
squares) with baryon number A = 5 at temperature T = 50MeV. 

Figure 13. The same as figure 12, but for baryon number A = 10. 

Figure 14. The same as figure 12, but for temperature T = 20MeV. 
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